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The 4-dimensional effective theory arising from an induced gravity action for a co-dimension greater 
than one brane consists of multiple galileon fields it 1 , I — 1, . . . , N, invariant under separate Galilean 
transformations for each scalar, and under an internal SO(N) symmetry. We study the viability of 
such models by examining spherically symmetric solutions. We find that for general, non-derivative 
couplings to matter invariant under the internal symmetry, such solutions exist and exhibit a Vain- 
shtein screening effect. By studying perturbations about such solutions, we find both an inevitable 
gradient instability and fluctuations propagating at superluminal speeds. These findings suggest 
that more general, derivative couplings to matter are required for the viability of SO(N) galileon 
theories. 



I. INTRODUCTION 

There has been much recent interest in theories of 
gravity arising from scenarios with extra spatial dimen- 
sions. Many examples of these are based on the Dvali- 
Gabadadze-Porrati (DGP) model [TJ [5] — a 4 + 1 dimen- 
sional theory with action consisting simply of separate 
Einstein-Hilbert terms in the bulk and on a codimension- 
1 brane, to which standard model particles are also con- 
fined. The model results in a 4D gravitational force 
law at sufficiently small scales, which transitions to a 
5D gravitational force law at a crossover length scale 
r c ~ Mpj /Mg, determined by the 5D and 4D gravita- 
tional couplings M 5 and Mp\ respectively. To yield in- 
teresting cosmological dynamics, this crossover scale is 
usually chosen to be of order the horizon size. 

Much of the phenomenology of the DGP model is cap- 
tured by its decoupling limit Mpi, M5 — > 00 with the 
strong-coupling scale A 5 ~ Mf /Mpi kept fixed [312]. In 
this limit, the difference between DGP gravity and gen- 
eral relativity is encoded in the behavior of a scalar degree 
of freedom, 7r. The dynamics of this scalar are invariant 
under internal galilean transformations ir — > n + c+b^x^, 
with c a constant and a constant vector. This sym- 
metry proves to be extremely restrictive, with a leading 
order self-interaction term which is a higher-derivative 
coupling cubic in tt, and yet which yields second or- 
der equations of motion. Higher order couplings with 
these properties were derived independently of the DGP 
model [SHE] and dubbed "galileons." See [5HTT] for cos- 
mological studies of galileon theories. 

It is natural to explore induced gravity models in co- 
dimension greater than one [18H26) . and recently multi- 
galilcon actions arising in the relevant 4-dimensional de- 
coupling limit have been derived [27 - 30 . The theories 
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studied in |30| are invariant under individual galilean 
transformations of the 7r fields, and also under an inter- 
nal SO(N) symmetry rotating the fields into one another, 
thus forbidding the existence of terms containing an odd 
number of n fields, in contrast to the co-dimension one 
DGP case. 

In this paper we explore the nature of spherically sym- 
metric solutions in theories with an SO(N) internal sym- 
metry among the galileon fields, and couplings to mat- 
ter that respect this symmetry. Spherical solutions for 
a more general bi-galileon action were discussed in (31j . 
for the specific case of a linear coupling ~ ttT to matter, 
where T is the trace of the matter energy momentum ten- 
sor. This form of coupling arises from decoupling limits of 
DGP-like theories, because tt arises through a conformal 
mixing with the graviton. However, while this coupling 
is therefore the natural form to consider in the case of a 
single galileon field, it breaks the new internal symmetry 
satisfied by multiple galileons (and breaks the galilean 
symmetry if the matter is dynamical) . We instead study 
general non-derivative couplings to matter fields which 
respect the SO(N) internal symmetry. 

At the background level, our solution can always be 
rotated to lie along a single field direction, say m , while 
the other field variables remain trivial, thus exhibiting 
spontaneous symmetry breaking. The solution exhibits 
Vainshtein screening [32] [33], characteristic of galileon 
theories: we find 7Tx ~ r sufficiently close to the source, 
whereas ~k\ ~ 1/r far away, with the crossover scale de- 
termined by a combination of the galileon self-interaction 
scale and the coupling to the source. However, when 
we turn to the stability of spherically symmetric solu- 
tions under small perturbations, we find that, sufficiently 
close to the source, perturbations in 7ri suffer from gra- 
dient instabilities along the angular directions. Morever, 
they propagate superluminally both along the radial and 
angular directions (in the regime that angular perturba- 
tions are stable). Perturbations in the remaining N — 1 
galileon fields are stable but propagate superluminally in 
the radial direction. 

The gradient instability and superluminal propaga- 
tion found here for the tti field are multi-field general- 
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izations of single galileon instabilities 5J. Our findings 
thus present significant hurdles for SO(N) galileon mod- 
els with non-derivative matter coupling. One of the main 
lessons to be drawn is that more general matter cou- 
plings, including derivative interactions, are necessary for 
the phenomenological viability of SO(N) mult i- galileon 
theories. For instance, the coupling ~ d^w 1 d v -KiT^ v nat- 
urally arises from brane- world constructions [301 33] and 
maintains both the galilean and the internal rotation 
symmetries. 



II. THE MODEL 

In co-dimension N, the 4-dimcnsional effective theory 
contains ./V fields tt 1 , I = 1---N, representing the N 
brane-bending modes of the full 4 + iV-dimensional the- 
ory. The extended symmetry of the vacuum lagrangian 
is 



Sir 1 = 



(1) 



where oj 1 ^, e 1 and uj 1 j are constant transformation pa- 
rameters. (See [30] for the geometric setup and origin 
of this symmetry). This transformation consists of a 
galilean invariance acting on each of the tt 1 fields, and an 
SO(N) rotation symmetry under which 717 transforms as 
a vector. The unique four dimensional lagrangian density 
respecting this is [23 130] 



r 
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-A [d^dv-ir 1 (d x d^njd x d l, 7r I - ^^tt/Dttj)] , 

(2) 

where A is a coupling with dimension [mass] -6 , contain- 
ing the strong interaction mass scale. The /, J indices 
are raised and lowered with 5jj. 

It remains to couple this theory to matter. The nat- 
ural coupling we might consider, the lowest dimension 
coupling that preserves the galilean and internal rotation 
symmetries, is ~ d^ir 1 d u iriT^ . This is the coupling that 
naturally arises from brane matter in the construction 
of [301 EH]- However, for static non-relativistic sources 
T^ u ~ P^o^O' an d since Oq-k = for static solutions 
there are no nontrivial spherically symmetric solutions 
with this coupling. 

Linear couplings Ainear ~ arise naturally from 
DGP-like setups, since the 7r's conformally mix with the 
graviton. These lead to spherical solutions [3T], but break 
the SO(N) internal symmetry. 

We therefore do not consider these couplings fur- 
ther, and instead concentrate on the most general non- 
derivative coupling that preserves the SO{N) symmetry. 
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where P is an arbitrary function of the invariant 7r 



III. SPHERICALLY SYMMETRIC SOLUTIONS 

Our focus is on the existence and viability of spher- 
ically symmetric solutions sourced by a delta function 
mass distribution 1 



T = -MS 3 (r 



(4) 



The equations of motion, including the coupling ([3|, are 
Dtt 1 - A [On 1 (d fl d u Trjd t *d' J ir J - On j Dttj) 

= AfP'(^ 2 (0))^(0)<S 3 (rO , (5) 

where P'(X) = dP/dX. Restricting to spherically sym- 
metric configurations 7r J (r), this reduces to 



1 d 

r 2 dr 

where 



[r 3 (y 1 + 2Xy I y 2 )] = MP' (tt 2 (0)) 7r / (0)J 3 (r) 
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r dr 



(6) 



(7) 



and y 2 = y I yi- Note that, due to the shift symmetry of 
the lagrangian, the equations of motion of galileon fields 
always take the form of a total derivative. Thus we can 
integrate once to obtain the equations of motion 



y 1 + ZXy'y 2 



M 
Aitr 3 



P' (7T 2 (0))7T'(0) 



(8) 



Dividing these equations by each other, we obtain the 
relations 



d-w 1 /dr _ 7r J (0) 
d7r J /dr ~~ 7r J (0) ' 

which, when integrated from the origin, gives 
7r x (r) _ tt j (0) 

7T J (r) ~ 7T J (0) ' 



(9) 



(10) 



The various components of the solution are therefore 
always proportional to each other. Thus, by a global 
SO(N) rotation, we can rotate the solution into one di- 
rection in field space, say the I — 1 direction, so that the 
solution takes the form ir 1 = ir and tt 1 = for / ^ 1. This 
model therefore exhibits a kind of spontaneous symme- 
try breaking of the internal SO{N) symmetry, since any 
non-trivial solution must pick a direction in field space. 
Equation (l8| now takes the form 



2Ay 3 = 



M 
4 7rr 3 



p> (^ 2 (o)Ho) 



(ii) 



1 Note that stable, non-trivial solutions without a source do not 
exist |35| . 
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As r ranges from zero to infinity, the left hand side is 
monotonic, and is positive or negative depending on the 
sign of P' (7r 2 (0)) 7r(0). For there to be a continuous so- 
lution for y as a function of r, the left hand side must be 
invertible when it is positive (negative). For a solution 
to exist, this requires (for non-trivial A) 



7T 



A > . 



(12) 



Thus y is also positive (negative), is monotonic with r, 
and ranges from zero to (negative) infinity as r ranges 
from infinity to zero. This in turn implies that d7r/dr 
does not cross zero, and hence 7r is monotonic. 



Equation ( 11 ) yields a solution for y, and hence d7r/dr, 
as a function of r and the parameters of the theory. In- 
tegrated from r = to infinity, this will give a relation 
between n(0) and the asymptotic value of the field 7r(oo). 
The asymptotic field value is essentially a modulus of 
the theory — it will be set by whatever cosmological ex- 
pectation value is present. It is a physically meaningful 
parameter as it affects the coupling to the source by de- 
termining 7r(0). 

Near the source, where the non-linear term dominates, 
the solution is linear in r, 
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FIG. 1: Schematic sketch of the solution for n(r) 



We expand the field in perturbations around the back- 
ground solution 7Tq, 



flr<r„(r) ~ tt(0) + 



M 
8^A 



P' (7r 2 (0))7r(0) 



1/3 



(13) 



whereas far from the source, where the linear term dom- 
inates, the solution goes like 1/r, 

M 1 

n r »rAr) ~ tt(oo) - ^P' (*») vr(0)- , (14) 

where the transition between these regimes occurs at the 
radius 



(AM 2 [P' (tt 2 (0))^(0)] ; 



1/6 



(15) 



Note that this crossover radius, and hence the distance 
at which non-linearities become important, depends on 
the modulus 7r(0). The equation of motion for ir(r) is 
readily solved numerically, and the solution obtained is 
plotted schematically in Fig. [T] 



IV. PERTURBATIONS: STABILITY AND 
SUBLUMINALITY 

While the existence of static, spherically-symmetric 
configurations is encouraging, there are, of course, other 
important checks that our solution must pass to be phys- 
ically viable. Specifically, following [5], we must study 
the stability of these spherically symmetric solutions and 
to determine the speed at which fluctuations propagate, 
since superluminal propagation can be an obstacle to 
finding an ultraviolet completion of the effective the- 
ory [55] , 



tt 1 = 7Tq + Sir 1 



(16) 



Away from the source, the linearized equations of motion 
for the perturbations are of the form 



-Kfr^S^+^dr (r 2 ^;(r)9 r (57r / )+^p(r)a^7r / = , 

(17) 

where the coefficients Kj(r), K'j(r) and Kf(r) depend 
on r through the background field tt 1 (r). We find 



K\ = 
Kl = 
K? = 



J_ d_ 
3r 2 dr 
1 + 6Ay 2 ; 

2r dr L 

1 d 
3r 2 dr 
1 + 2\y 2 

1 d 
2r dr 



[r 3 (1 + 18Ay 2 )] 
, 2 ; 

r 2 (l + 6Xy 2 )] ; 

[r 3 (1 + 6Xy 2 )] ; 
? ■ 

r 2 (l + 2Xy 2 )] . 



(18) 



Applying the implicit function theorem to the function 
F(y,r) =y + 2Xy 3 - ^P' (tt 2 (0)) tt(0) = 0, we have 



dy 
dr 



d r F 3 y + 2Xy 3 



d y F r 1 + 6Ay 2 



(19) 
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This allows us to eliminate dy/dr from (18 1: 



K\ = 


(1 


-6Ay 2 ) 2 


1 


+ 6Xy 2 


K{ = 


1 + 


GXy 2 ; 


= 


1 - 
14 


6At/ 2 
6Xy 2 ' 


Kj^ = 


14 


12A 2 y 4 


1 


+ 6Xy 2 


K r j^ = 


1 + 


2Ay 2 ; 


K% x = 


14 
14 


2Xy 2 
6Xy 2 



(20) 



Stability of the spherically symmetric background so- 
lutions against small perturbations requires K > for 
all K's. The I 1 directions in field space are stable, 
but the 7T 1 direction exhibits a gradient instability suffi- 
ciently close to the source along the angular directions. 
In other words, < near the source. Therefore, lo- 
calized perturbations can be found near the source that 
lower the energy of the solution through their gradients. 
This instability plagues very short-wavelength fluctua- 
tions, right down to the UV cutoff, so decay rates are 
dominated by the shortest distances in the theory and 
cannot be reliably computed within the effective theory. 



Equations (20 1 also allow us to compute the speeds 
of propagation of our small perturbations, in both the 
radial and angular directions. These are given by 



[c 2 )l 



K{ _ /14-6A?/ 2 
~K\ _ Vl-6Ay 2 



1 



1 

-6Ay 2 ' 

1 + 2Ay 2 ) (l 4- 6Ay 2 ) 
1 + 12A V 

-2Xy 2 
12AV ' 



1 



(21) 



Note that (c 2 )^ > 1, and hence these perturbations al- 
ways propagate supcrluminally. The same is true of 
(c 2 )f , in regions where these perturbations are stable. 
The speed (c 2 )J_^ 1 is always superluminal, and (c 2 )p_^ 1 is 
always subluminal. Whether superluminal propagation 
of signals is problematic for a low-energy effective theory 
is still an arguable issue, but it seems that at the least it 
may preclude the possibility of embedding the theory in 
a local, Lorcntz-invariant UV completion [36 . 



A. Other constraints 

It is interesting to note in passing that if a mechanism 
exists to tame the instabilities we have identified, then 
precision tests of gravity within the solar system already 



place useful constraints on multi-galileon theories. The 
galileon is screened at radii below the Vainshtein radius 
r», given by equation (fl5|), restoring the behavior of gen- 



eral relativity. Requiring the solar system to be screened 
to r ~ 10 16 m thus yields a constraint on A and 7r(0). 
However, lunar laser ranging data constrain the depar- 
ture from the gravitational potential predicted by GR to 
satisfy f < 2.4 x lO" 11 (at radius r = 3.84 x 10 10 cm), 
and we may translate this into a constraint on a different 
combination of A and 7r(0) 

For example, consider the choice of P(X)Mpi ~ 
\J w i tvi, giving a linear coupling between the radial tt field 
and matter. In the interesting case when the constraints 
are saturated, and detection of an effect is therefore im- 
minent, the relevant constraint simply becomes 



1 



A 1 / 6 



< 10~ 9 eV . 



(22) 



Note that this is an extremely low cutoff for the effective 
theory, as is also found in the DGP model. 



V. DISCUSSION 

We have derived spherically-symmetric solutions in an 
SO(N) multi-galileon theory with general, non-derivative 
couplings to matter. These solutions exhibit a Vain- 
shtein screening effect, characteristic of galileon models. 
However, a study of the behavior of fluctuations around 
these solutions shows that one of the fields has imagi- 
nary sound speed along the angular directions, signaling 
an instability to anisotropic modes of arbitrarily short 
wavelength. Moreover fluctuations inevitably propagate 
superluminally. 

These results raise serious concerns about the phe- 
nomenological viability of SO(N) multi-galileon theories. 
(Of course, this does not preclude their effectiveness in 
early universe physics [151 137j , for instance during infla- 
tion, as long as they become massive or decouple before 
the present epoch.) A key input in our analysis is the re- 
striction to non-derivative coupling to matter. The main 
lesson to be drawn is that more general, derivative cou- 
plings are necessary. For instance, the lowest-dimensional 
coupling invariant under the galilean and internal rota- 
tion symmetries is ~ d^n 1 d v iiiT^ v . This coupling in 
fact naturally arises in the higher-codimension brane pic- 
ture [34]. As mentioned earlier, the galileon fields are 
oblivious to static, spherically-symmetric sources in this 
case; thus exhibiting a screening mechanism. However, 
they will be excited by orbital motion, and we leave a 
study of the phenomenological implication of this cou- 
pling to future work. 

Our analysis also highlights a distinct advantage to ex- 
plicitly breaking the symmetry ([T]), for example through 
the introduction of a sequence of regulating branes 
of different co-dimensions, as in the cascading gravity 
case [HlllllMI- The explicit breaking of SO(N) sym- 
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metry allows for more general terms in the action, which 
can lead to a healthier phenomenology [31 j - 

Finally, should a creative cure for our instabilities be 
found, then we have demonstrated that precision so- 
lar system tests of gravity set interesting constraints on 
multi-galileon theories. 
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